Universal model for exoergic bimolecular reactions and inelastic processes 
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From a rigorous multichannel quantum-defect formulation of bimolecular processes, we derive a 
fully quantal and analytic model for the total rate of exoergic bimolecular reactions and/or inelastic 
processes that is applicable over a wide range of temperatures including the ultracold regime. The 
theory establishes a connection between the ultracold chemistry and the regular chemistry by show- 
ing that the same theory that gives the quantum threshold behavior agrees with the classical Gorin 
model at higher temperatures. In between, it predicts that the rates for identical bosonic molecules 
and distinguishable molecules would first decrease with temperature outside of the Wigner threshold 
region, before rising after a minimum is reached. 



PACS numbers: 34.10.+x,34.50.Lf,34.50.Cx,03.65.Nk 

The recent experiment by the JILA group [l( repre- 
sents a milestone in studies of chemical reactions. For the 
first time, reactions are studied in a temperature regime 
where the quantum nature of the relative motion of the 
reactants becomes unequivocally important, as reflected 
in the quantum threshold behavior and in the importance 
of quantum statistics. More importantly, the experiment 
strongly suggests that bimolecular reactions in the ultra- 
cold regime follow universal behaviors determined by the 
long-range interaction, as spelled out in more detail in 
related theoretical works by Julienne and Idziaszek 0, Q 
and by Quemener and Bohn 

The significance of such experiments goes beyond ex- 
ploring chemical reactions in a new temperature regime 
with many unique characteristics, such as controllabil- 
ity via moderate external fields 4-6]. By forcing a new 
perspective on the quantum theory of reactions, as de- 
manded by their interpretation, they have potential to 
improve our understanding of reactions and inelastic pro- 
cesses at all temperatures. This paper is a illustra- 
tion of such an outcome. From a rigorous multichannel 
quantum-defect formulation of bimolecular processes, we 
derive here a fully quantal model for the total rate of 
exoergic bimolecular reactions and/or inelastic processes 
that is applicable over a wide range of temperatures. 
The theory establishes a connection between the ultra- 
cold chemistry and the regular chemistry by showing that 
the same theory that gives the quantum threshold be- 
havior [lrli| agrees with the classical Gorin model 0, IH 
at higher temperatures. In between, it shows that the 
rates for identical bosonic molecules and distinguishable 
molecules would first decrease with temperature outside 
of the Wigner threshold region, before rising after a mini- 
mum is reached. The theory further illustrates explicitly 
how the quantum effects, including effects of quantum 
statistics, gradually diminish at higher temperatures, and 
establishes the van der Waals temperature scale as the 
one that separates the quantum and the semiclassical be- 
haviors of reactions. The same formalism is applicable to 
ion-molecule reactions where our quantum model, with 



details to be presented in a separated publication, would 
approach the classical Langevin model [1, 0] at high tem- 
peratures. 

Consider the collision of two distinguishable molecules 
A and B in the absence of any external fields. The cross 
section for a transition from an entrance channel i to an 
exit channel / can be written in terms of the S matrix 
as 10] 
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(2F Al + l)(2F Bl + l)kf 
x £ (2F t + l)|4f 

F t ,F„h,{q f } 



(1) 



Here e = E — = h kf/2[i is the energy relative to 
the entrance channel i, with \x being the reduced mass. 
Fa and Fb are the total (internal) angular momenta of 
molecules A and B, respectively. F = Fa + Fb is the to- 
tal angular momentum excluding I, which is the relative 
angular momentum between A and B. F t is the total an- 
gular momentum of the system, which is conserved in the 
absence of external fields. {(?/} represents the quantum 
numbers, excluding Ft, that are required to characterize 
an exit channel /. 

The exit channels can be classified into elastic chan- 
nels, labeled by {e}, inelastic channels, labeled by {u}, 
and reactive channels, labels by {r}. From the unitar- 
ity of the S matrix [10(, the total cross section, cr ur = 
r } a fii f° r ^ ne combination of all inelastic and re- 
active processes, can be written as 

CTur(e) 



{2F Al + l)(2F Bi + l)kf 
x £ (2F t + l)(l-X;iS<f | a ). (2) 



F t ,Fi,U 



F e ,l e 



The implication is that such a total cross section is com- 
pletely determined by the S matrix elements for elastic 
channels only. The corresponding rate constant at tern- 
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perature T is given in terms of <7 ur by 



K(T) = 



TTfJ, 



1/2 



(k B ry 



ecr ur (e) exp(— e/fcsT)(ie . 

(3) 

where ks is the Boltzmann constant. 

Considerably further understanding of bimolecular 
processes can be achieved through a multichannel 
quantum-defect theory (MQDT) (see Ref. [11( and ref- 
erences therein), especially through an 5* matrix formu- 
lation in terms of quantum reflection and transmission 
amplitudes associated with the long-range potential [l2| ■ 
The theory, which is a multichannel generalization of the 
S matrix formulation of Ref. [l2| , gives 

S ( *> =-(-!)' Hf+t^S^-r^S^-H^] , (4) 



■(-1)' <rg>+tMS° e 
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Here (—1)' is a diagonal matrix with elements (— 
for channel j. The rif and too*' are diagonal matri- 
ces for the open channels with elements rj (e s j) and 

t\ m \e S j) representing the (complex) quantum reflection 
and the quantum transmission amplitudes, respectively, 
for molecules going outside-in (approaching each other) 
[l2j ]. They are universal functions of scaled energies, 
e s j = (E — Ej) j SEj, that are uniquely determined by the 
exponent, ay, of the long-range interaction, —Caj/Rj 3 , 
in channel j, and the lj. Such long-range interactions 
have length scales /3 a j = (^f^jC a j/h' 2 ) 1 ^ aj ~ 2 ^ and cor- 
responding energy scales SEj = (ft 2 /2^j)(l//3^ ■), associ- 



ated with them. The 



and too ^ are similar, except 



that their elements are amplitudes for molecules going 
inside-out (moving away from each other) [ljj]. The S^ s 
is an effective short-range S matrix [JJj], after the elim- 
ination of the closed channels [H| . It has the physical 
meaning of being an effective reflection amplitudes by 
the inner potential. 

Equation Q for the S matrix has a clear physical in- 
terpretation as discussed in Ref. [12J. In particular, the 
rn-th term in its expansion, Eq. ([5]), corresponds to the 
contribution from a path in which the fragments are re- 
flected m+1 times by the inner potential. Further simpli- 
fication can be achieved by recognizing that rf°^ (e s j) ~ 
for e S j >• SEj [HI- Dividing all open channels into elas- 
tic and near-degenerate channels with e S j <~ SEj, and 
other channels with e S j 3> SEj, we have 



r {oi) S 
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4f} , (6) 



channels for which the quantum reflection amplitude r\ l °^ 
differs substantially from zero. 

A number of different theories and models, both exact 
and approximate, can be derived from Eq. @, and ei- 
ther Eq. (UJ) or ©. The universal model to be presented 
here, which we call the quantum Langevin (QL) model, 
results from the assumption of no reflection by the inner 
potential, namely, 



Sl s « 



(7) 



It is a rigorous mathematical representation of the 
Langevin assumption 0, Q in a quantum theory. In plain 
language, it assumes that whenever two molecules come 
sufficiently close to each other, so many "bad" things can 
happen that they can never get out of it in their initial 
configurations. It can be expected to be a good approx- 
imation whenever there are a large number of open exit 
channels that are strongly coupled to the entrance chan- 
nel in the inner region. For it to be satisfied in the limit 
of zero energy, the reactions and inelastic processes under 
consideration have to be at least exoergic. 

Under the Langevin assumption, Eq. ([6]) gives 
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(Ft) 



-(-1) 



(8) 



It implies that the elastic S matrix elements in the QL 
model, and therefore the total cross section and the cor- 
responding total rate for inelastic and reactive processes, 
are all described by universal functions that are uniquely 
determined by the long-range interaction in the entrance 
channel. Substituting Eq. © into Eq. @ and subse- 
quently into Eq. ([3]) , the total rate of inelastic and reac- 
tive processes in the QL model can be written as 



K{T) =a K KS a \T„) 



(9) 



Here sk is the rate scale corresponding to the long-range, 
—C a /R°, interaction in the entrance channel. 



(10) 



in which h/ fi/3 a is the velocity scale corresponding to the 
length scale (3 a . /C (q) (T s ) is a universal function of the 
scaled temperature, T s = T/(s£://cb), that is uniquely 
determined by the exponent a. Specifically, 



JC {a) {T s 



dxx 1 ' 2 e- x yV [a \T s x) , (11) 



where W {a) {t s ) is a scaled total rate before thermal av- 
eraging. It depends on energy only through the scaled 
energy e s = e/sE, and has contributions from all partial 
waves: 



where S£ s is a submatrix of the effective short-range S^g 
that includes only the elastic and other near-degenerate 



W( Q )(£ s ) = ^W ; (a) ( £s ) 



(12) 
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Here is a scaled partial rate given by 



W i (Q) (e s ) = (2/ + l)7; c(a) (e s )/ey 2 



(13) 



in which T^ a \t s ) = \t\ m, (e s )\ 2 is the quantum transmis- 
sion probability through the long-range potential at the 
scaled energy e s and for partial wave I 12j. 

This QL model for reactions and inelastic processes 
is applicable to both neutral-neutral systems, for which 
a = 6 corresponding to the van der Waals potential, and 
charge-neutral systems, for which a = 4 corresponding to 
the polarization potential. We focus here on the neutral- 
neutral case to make connection with existing theories 
and experiments [lJ-Q] - The results for charge- neutral 
systems will be presented elsewhere. 

For a = 6, the quantum transmission probability 
through the long-range potential, 7j (e a )j which is the 
only quantity required to determine the universal rate 
functions in the QL model, can be found analytically by 
substituting Eqs. (Al)-(A4) of Ref. into the Eq. (52) 
of the same reference. The result is 



Aoi) 



7T (6) ( £s ) 



2M ts i \cos(irv) — cos(3nv)] 
1 - 2M eJ cos(3tw) + Af 2 \ 



(14) 



Here v is the characteristic exponent for —1/R 6 potential 
[H, and 



M eJ (v) = |A| : 



r(i- v ) 



r(i 



l/Q-V 



r(i + VQ - v) 



|r(i - 



r(i + v + v 

~c«.i(-")~ 



(15) 



(16) 



where A = e s /16, v a = (21 + l)/4, and 

CO 

c e .i(u) = llQ(y + j), 

3=0 

in which Q(v) is given by a continued fraction: 



(17) 

The resulting universal rate function, KS°^(T S ), applica- 
ble to neutral- neutral distinguishable molecules, is illus- 
trated in Figure [T] Similar results can be obtained for 
neutral- neutral interactions of identical molecules, fol- 
lowing considerations similar to those of Ref. [IH • They 
are given generally by a combination of two universal rate 
functions defined by 



JC s{a) (T s 



dxx 1/2 e- x W s(a) (T s x) , (18) 



where 



i-even 



VV/ a) (e s ) 



(19) 
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FIG. 1. (Color online) The universal rate functions K, s(a) (T 3 ), 
tC A{a) (T s ), and JC {a) (T s ) for a = 6, corresponding to -1/R 6 
type of interaction in the entrance channel. Here /Cg refers 
to the prediction of the classical Gorin model, as given by 
Eq. (1251) . All results include summations over all relevant 
partial waves. 



and 



9 1"°° 

)CM<*)(T S ) = — dxx 1/2 e- x W A{a) {T s x) , (20) 
Jo 



where 



W A ^(e s )=2 W z (Q) (e s ) 



(21) 



/=odd 



For example, in terms of K. s{a) and JC Ma \ K S (T) = 
skK- S ^ gives the rate for identical bosonic molecules in 
the same internal (M) state, and K A (T) = s K JC A{ - a) 
gives the rate for identical fermionic molecules in the 
same internal state. The three rate functions are related 
by /C( Q > = (/C S ( Q ) + K A ^)/2, and are all illustrated in 
Figured] for a = 6. 

At ultracold temperatures such that T s « 1, a QDT 
expansion [HI of (e s ) gives 



JC S ^(T S ) = 8a sl=0 



4a. 



M T y 2 + 3 4 =0 T s + O(T s 3 / 2 ) 



where a si=0 = 2vr/[r(l/4)] 
mean scattering length for I = 15| . 

16a s z = i 



(22) 

0.4779888 is the scaled 



JC A ^{T S ) = 36a sl=1 T s 



1 



T s 3/2 + 0(T 2 ) 



(23) 

where a sl=1 = [r(l/4)] 2 /367r w 0.1162277 is the scaled 
mean scattering length for I — 1 [15j . and 



/C (6) (T S 



._ (4a s ;=o) 2 rr i/2 
4a s / =0 7= — T s 



+ (l2a 3 sl=0 + 18o.,=i) T s + 0(T^ 2 ) .(24) 
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At high temperatures as characterized by T s 3> 1, it is 
straightforward to show, from the semiclassical limit of 
the transmission probabilities [l2j], that 



K 5(6) (T S ) « /C A(6) (T S ) a /C (6) (T S ) 



2 4 / 3 r(2/3) 



T l/6 



(25) 

in agreement with the classical Gorin model [7J, |8[ . All 
scaled results can be put on absolute scales using a single 
parameter, the C§ coefficient for the entrance channel, 
from which both the temperature scale Ss/fcs and the 



rate scale sk can be determined 

In the Wigner threshold region, in which the rates are 
accurately characterized by the first terms of Eqs. (|2"2")l - 
([24]) . our results are consistent with those of Julienne 
and Idziaszek 2, [§]. Outside of this region, both IC^ 
and /C s < 6 ) are predicted to first decrease with tempera- 
ture, a behavior that deviates strongly from the predic- 
tion of the classical Gorin model. Specifically, K.^ is 



predicted to reach a minimum value of K, 



(6) 



1.587 at 



T 



(6) 



0.1154, for a drop of about 17% from its value at 
zero temperature. The K, si ^> is predicted to reach a min- 
imum value of « 1.908 at « 1.114, for a drop 
of about 50% from its value at zero temperature. For 
the JILA experiment [l[, T^ in translates, using the Cq 
coefficients of Kotochigova 



16, 171, to T 



(6) 



2.58 fiK 

for 40 K 87 Rb+ 40 K 87 Rb in different internal states, and 
to T^) n w 11.9 [iK for 40 K+ 40 K 87 Rb. It is worth not- 
ing that an experimental measurement of either 



or 



^min" 1 wou ld constitute a measurement of the Cq coef- 
ficient, a fact that can be valuable especially for more 
complex molecules for which theoretical calculations of 
C 6 0, [H become increasingly difficult and unreliable. 
At higher temperatures, our results show how the quan- 
tum effects, including that of quantum statistics, gradu- 
ally diminish, and all rates approach that of the classical 
Gorin model @, HI ■ As illustrated in Fig. (TJ such a transi- 
tion from quantum to semiclassical behavior occurs over 
a range of the van der Waals temperature scale s^/fca. 

The QL model gives the total rate that includes both 
reactive and inelastic processes. For experiments with 
only reactive channels open [H, it give the total rate of 
reactions. For experiments with no open reactive chan- 
nels, it gives the total rate of inelastic processes. In all 
case, the requirement for its validity is that there are 
many open channels that are strongly coupled to the en- 
trance channel by the short-range interactions. Of inter- 
est in the context of cold-atom physics, the QL model 
serves to unify theories of ultracold chemistry Il7 | 
with theories for atom-atom [lj| , atom-molecule |2of and 
molecule-molecule inelastic processes. For example, for 
a vibrational highly excited molecule, except Feshbach 
molecules with with very small binding energies 2l| , the 
theory predicts that its collisional lifetime is approxi- 
mately independent of its initial state, and that the rates 



for atom-molecule and molecule-molecule inelastic pro- 
cesses are related. More specifically, for Cs2 in a highly 
excited rovibrational state [13] , it predicts that the Cs 2 - 
Cs2 inelastic rate should have a minimum around 6.31 
/iK (assuming that they are prepared in the same state) , 
and Cs2-Cs inelastic rate has a minimum around 1.70 
/iK. More detailed discussion of such applications will be 
presented elsewhere. 

In conclusion, we have presented a universal model of 
exoergic bimolecular reactions and/or inelastic processes 
that is applicable over a wide range of temperatures, il- 
lustrating the evolution from quantum behavior to semi- 
classical behavior. It is an important baseline model in 
which rates for different systems differ from each other 
only in scaling, and has an intriguing and useful property 
of being more accurate for more complex systems. Simple 
analytic formulas, to higher orders than those of Julienne 
and Idziaszek [2, 3], and applicable over a substantially 
wider range of temperatures, are also presented. Equally 
important, we believe, is that the underlying MQDT for- 
mulation, such as Eq. (j6|), lays a solid foundation for new 
types of theories of reactions and inelastic processes, ei- 
ther rigorous or approximate, that goes further beyond 
the QL model. 
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